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Abstract. A Wigner crystal structure of the electronic ground state is induced by strong Coulomb interac- 
tions at low temperature in clean or disordered two-dimensional (2d) samples. For fermions on a mesoscopic 
disordered 2d lattice, being closed to a torus, we study the persistent current in the regime of strong in- 
teraction at zero temperature. We perform a perturbation expansion starting from the Wigner crystal 
limit which yields power laws for the dependence of the persistent current on the interaction strength. The 
sign of the persistent current in the strong interaction limit is independent of the disorder realization and 
strength. It depends only on the electro-statically determined configuration of the particles in the Wigner 
crystal. 



PACS. 73.23.Ra Persistent currents 
localized states 



71.10.Fd Lattice fermion models - 71. 23. An Theories and models; 



1 Introduction 

The interplay of disorder and interaction in mesoscopic 
samples has attracted considerable interest in recent years 
M. A prominent experimental finding in this field is the 
insulator-metal transition in 2d systems of high mobil- 
ity g|, occurring when the carrier density is increased. 
This transition can neither be explained by disorder ef- 
fects alone nor by interaction effects in clean samples and 
is the subject of still increasing experimental and theoret- 
ical activities (for a review see || and references therein). 

One possible mechanism invokes both, interaction and 
disorder and associates the insulator-metal transition 
with the melting of a pinned Wigner crystal when the 
carrier density increases and therewith r s (the interaction 
energy in units of the kinetic energy) decreases, allowing 
for metallic behavior at intermediate r s . Only at a higher 
density, when r s is small, the interaction becomes negli- 
gible and the now dominating disorder-induced Anderson 
localization leads back to insulating behavior. Such a sce- 
nario is supported by experiments B and numerical inves- 
tigations of few interacting particles in disordered lattice 
models |,|,[7|. 

Another important experimental result in mesoscopic 
physics whose explanation necessitates to invoke disor- 
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der and interactions simultaneously is the value of the 
persistent current in diffusive rings |8|. These persistent 
currents are much larger than the theoretical prediction 
for non-interacting electrons in disordered rings Q . While 
the electron-electron interaction seems to play an essential 
role, the disorder in the sample is also important: Inter- 
actions cannot affect the persistent current in clean rota- 
tionally invariant Id rings ]To|,pT[, and the non-interacting 
result is consistent with the experimental one for a clean 
semiconductor ring in the ballistic regime fl2]| . 

This has generated a large theoretical activity, dealing 
with the combined effect of interactions and disorder on 
the enhancement of persistent currents in mesoscopic rings 
(for an overview see e.g. Jj],[l3|,[l4[ and references therein). 
Even though different theoretical approaches suggest an 
increase of the persistent current in disordered samples 
due to repulsive Coulomb interactions, a quantitative un- 
derstanding of the experiments is still lacking. 

Within continuous models of fermions in Id with dis- 
order, repulsive interactions are found to enhance the per- 
sistent currents, without jllj] and with fl5|| spin. On mod- 
erately disordered lattices, repulsive interactions are how- 
ever found to decrease the persistent current for spinless 
fermions in Id p6[ |l7|,a q|. It was concluded from analyt- 
ical considerations fll 9f , using renormalization group ar- 
guments, that Id lattice models exhibit an interaction- 
induced enhancement of the persistent current only when 
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the spin degree of freedom is taken into account. Never- 
theless, a numerical study of small systems ||(J revealed 
that repulsive interactions can slightly enhance the aver- 
aged persistent current even for spinless fermions on a Id 
lattice, provided the disorder is very strong. 

A recent numerical investigation treating the persis- 
tent current for strongly disordered individual Id chains 
pl|,p2| leads to the conclusion that, while Anderson local- 
ization is dominating the non-interacting case, the persis- 
tent current can be strongly enhanced by repulsive local 
interactions at sample dependent intermediate values of 
the interaction strength. This derealization is accompa- 
nied by a reorganization of the ground state structure. 
At half filling, strong interactions induce a regular charge 
density in the Mott-insulator regime and decrease the per- 
sistent current. 

This suppression of the persistent current in the limit 
of strong repulsive interactions is not limited to Id models 
of spinless fermions at half-filling with local interactions. It 
occurs with long-range interactions 1 16 at arbitrary filling 
and also in 2d lattices ||, when strong interaction (large 
r s ) leads to the formation of a Wigner crystal pinned by 
the disorder in the ground state. 

In the case of the Hubbard model at half filling (2^] and 
for spinless fermions in Id chains without |24| and with 
disorder p5[ , this suppression of the persistent current 
at strong interaction has been understood quantitatively 
from a perturbation theory starting at the Mott insula- 
tor limit. In contrast, numerical Hartree-Fock approaches 
which allow to treat larger systems than the exact diago- 
nalization used in Refs. |q,^q] and which have been used 
for weakly interacting particles in Id and 2d [Q, are, at 
least in Id, not able to quantitatively describe the persis- 
tent current in the limit of strong interaction [^7j . 

The sign of the persistent current in the case of spinless 
particles in strictly Id rings is independent on the disor- 
der realization and the interaction strength. According to 
a well-known theorem by Leggett |28| , the sign is given 
solely by the parity of the particle number N (paramag- 
netic for N even and diamagnetic for N odd) . This general 
rule is confirmed by an explicit calculation for a Luttinger 
liquid ring without disorder [ p9[ . 

On the other hand, the sign of the persistent current 
for particles with spin in disordered Id rings or spinless 
fermions in disordered 2d systems differs from sample to 
sample and no general sign rule exists. Only in some spe- 
cial cases, like non-interacting particles with spin in clean 
Id rings |30fl , and for some interacting situations using the 
Hubbard model |3l],|32|], the sign of the persistent current 
has been determined. The situation is less clear in the pres- 
ence of long-range interactions and disorder we address in 
this work. For electrons with spin in Id rings having a 
particular disorder consisting of only one barrier, a ten- 
dency towards diamagnetic responses, independent of the 
particle number, was found at strong repulsive interaction 



For spinless fermions in strongly disordered 2d lat- 
tice models, it has been noticed in numerical studies that 
the sign of the persistent current becomes realization- 



independent in the limit of strong interaction |5|,|6|,|34|,[35||. 
Detailed studies of the local current show that the sup- 
pression of its transverse component by the interactions is 
much stronger than the decrease of the longitudinal cur- 
rent. On such a lattice, closed to a torus, the structure of 
the ground state at strong interaction is a Wigner crystal 
pinned by the disorder ^,g,^5). While the system exhibits 
Anderson localization at weak interaction, the regime of 
intermediate interaction shows indications of a new type 
of correlated metal |3(| . 

For spinless fermions in 2d lattice models without dis- 
order, the amplitude of the persistent current has recently 
been studied analytically and numerically at strong inter- 
action |5?J . When r s is large, the hopping matrix elements 
between neighboring lattice sites being much smaller than 
the interaction strength, the behavior can be understood 
from a perturbation theory in terms of the hopping matrix 
elements. 

In this paper we report a study of the persistent 
current in disordered 2d lattice models at very strong 
Coulomb interaction, using a perturbation theory expan- 
sion around the pinned Wigner crystal. While the under- 
standing of this regime cannot directly explain the high 
persistent currents observed experimentally in diffusive 
metal rings, it may be relevant for the insulating side at 
low carrier density of the insulator-metal transition in 2d. 
We show how the sign of the persistent current at strong 
interaction follows systematically from the structure of the 
Wigner crystal and find simple rules for this sign. For the 
absolute values, power laws similar to the ones found in 
[ p7| are obtained also in the disordered case. We shall show 
how the presence of disorder and spin can influence the 
prefactors and the exponents of these power laws. 

In the following section we introduce the model for 
interacting fermions on a disordered lattice and the quan- 
tities used to characterize its properties. The perturbation 
theory is developed in section [| and applied in section ^ 
to the persistent current in longitudinal and transverse 
direction before we conclude the paper. 



2 Model 

2.1 Hamiltonian 

We consider N fermions on a disordered square lattice 
with Coulomb interaction. The corresponding Hamilto- 
nian reads 

H = H K + H D + H 1 . (1) 



The kinetic energy term is 



(2) 



with the hopping matrix element t = 1 setting the energy 
scale. We concentrate on rectangular 2d lattice structures 
with L x x L y sites i. The fermionic on-site operators Ci i(T = 
c {xi,yi),a destroy a particle with spin a located at r-j = 
(xi, yi), where the position coordinates Xi € {1, 2, . . . , L x } 
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and y l £ {1, 2, ... , L y } are measured in units of the lattice 3 Theoretical approach 

spacing a. 

The disordered potential contribution is 



(3) 



where W is the disorder strength, with the independent 
random variables Vi , drawn from a box distribution within 
the interval [—1/2; +1/2]. The occupation number oper- 
ators are as usual given by hi^ = c^c^ a . The Coulomb 
interaction is described by the term 



H 



b l,<7 'H' .<T 



U J2 



(4) 



the interaction strength being parametrized by U. The 
sum in H-& runs over all pairs of sites which are next 
neighbors on the lattice < >, while the interaction 
term is composed of a sum over all pairs of different sites. 



With these definitions, one finds r s = 



with the 



2t\ 

electronic density or filling factor v — N/L x L y = l/b 2 , 
b being the average distance between particles. An addi- 
tional term takes into account double occupancy of a site 
by two particles of different spin, d (< a) being a measure 
for the size of the on-site orbitals. 

In order to study the persistent current, we close the 2d 
lattice first to a cylinder by imposing generalized periodic 
boundary conditions 



(5) 



in ^-direction. For <f> x — 0, this is equivalent to usual pe- 
riodic boundary conditions. Finite <f> x accounts for a mag- 
netic flux <P — (f> x <Po/2ir threading the ring, 4>q being the 
flux quantum. We choose the units such that $o/27r = 1. 
In order to reduce finite size effects, we use periodic bound- 
ary conditions in y-direction, and thus obtain a t orus 
topology with the fluxes {4> x ,4> y ) — (4> x ,0). In section (4.3) 
we will use the dependence on the transverse flux <p y to 
study also the transverse current. 



3.1 Wigner crystal at strong interaction 

In the non- interacting limit, disorder leads to Anderson 
localization of the one-particle states and the problem can 
be treated by a perturbative expansion around the on-site 
localized states in terms of the hopping matrix elements t 
[p9| . Hopping to distant sites costs disorder energy of the 
order of W such that a series expansion in t/W results. 

In the many-body case, strong repulsive interaction U 
leads to Wigner crystallization of the ground state with 
on-site localized charges in the electro-statically most fa- 
vorable position. One can use a similar perturbative for- 
malism in terms of the hopping i, but now, the essential 
cost in energy caused by deplacing one of the many parti- 
cles is given by the increase of the interaction energy such 
that one obtains a systematic expansion in terms of t/U . 

We decompose the Hamiltonian (13) as 



H = Hq + Hi, 



(7) 



with an unperturbed part containing disorder and inter- 
action 

H = H D +Hj, (8) 

and the perturbation given by the hopping terms of Hk- 
Hq is composed of terms containing only occupation 
number operators in the one-particle on-site basis. There- 
fore, its A^-particle eigenstates \ip a ) are Slater determi- 
nants built from N different one-particle functions and 
are completely characterized by the occupation numbers 
n%,a{ot) G {0, 1} of the one-particle states on site i with 
spin er, fulfilling the condition N = J^. rij i(T (a). There- 
with, the many-body eigenstates of Hq can be written in 
the form 



r(ot) 



10) 



(9) 



(|0) is the vacuum state), and the corresponding eigenen- 
ergies are given by E a = + E l a with 



E a = W ^2 v i n i,"( a ) and 



(10) 



2.2 Persistent current 

The magnetic flux threading the ring can drive a persistent 
current through the system. At zero temperature, it is 
given by 

9E 



(6) 



where E is the many-particle ground state energy. Thus, 
the persistent current at T — is a measure of the depen- 
dence of the ground state energy on the magnetic flux. 
Since the latter can be expressed in the form of a bound- 
ary condition, it is at the same time a measure of the 
ground state sensitivity to the boundary conditions and 
can be related to the conductance of the sample [B8| . 



n it<7 {a)n v<cr i{a) 



(11) 



The ground state |-0o) of this Coulomb glass problem is 
given by purely classical considerations, minimizing dis- 
order and interaction energy. Its charge configuration de- 
pends in general on the specific disorder realization of the 
sample. At strong enough interaction, when the disorder 
effects are dominated by the interaction, the structure of 
\ipo) is the Wigner crystal of minimal interaction energy, 
independent of the disorder realization. The rigid array of 
charges can be translated as a whole through the system 
without changing the interaction energy. Nevertheless, the 
contribution of the disordered potential to the energy de- 
pends on the realization and pins the Wigner crystal in a 
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realization dependent position. It is important to realize 
that the structure of the Wigner crystal is entirely given 
by the lattice geometry and the Coulomb interaction of 
the N particles. 

In contrast to Hq, the perturbing part £/k of the 
Hamiltonian depends on the magnetic flux through the 
ring since the latter appears in the boundary condition 
in a;-direction and therefore influences some hopping ma- 
trix elements. Writing all the hopping terms explicitly, one 
obtains 

(L* L v 
E E C t,lOi* < W-l);cr 
x— 1 y— 2 



' E E C tx,y);<r C (x-l,y);a + £ ^ E C (*,l);<: 



-(x,L y );o 



x=2y=l 



-e'^Vct s c (T , +H.C. 



(12) 



3.2 Perturbation expansion for spinless fermions 

For electrons with spin, the unperturbed ground state is 
2 JV -fold degenerate since all the spin configurations yield 
the same energy when hopping is suppressed. Further de- 
generacies appear in the case of clean systems when W = 
allows for translational symmetry. 

For simplicity we first treat the case of completely 
spin polarized systems (all spins up, equivalent to spin- 
less fermions) with disorder where the ground state of Ho 
is not degenerate and the expansion in Hk is straightfor- 
ward. 

Using standard perturbation theory, the correction to 
the ground state energy in n th order is given by 



E 



(n) 



E 



(13) 



(ipn\H K \i/j ai ) (ip ai \H K \t/j a2 ) . . . (V>q„_ 1 |i?K|V'o) 
(Eq — E ai )(E — E a2 ) . . . (Eo — E otn _ 1 ) 

with the sums running over all the eigenstates of H except 
the ground state itself. 

The numerator of the contributions to the sum of equa- 
tion ( |l3| ) contains matrix elements (ip ak I-HkIVw+i) of the 
perturbing Hamiltonian. Since consists only of one- 
particle hopping terms, non-zero matrix elements can arise 
only if the two states \tpa k ) and \ip ak+1 ) differ by nothing 
else than a single hop of one of the particles. From (O) 
one sees that a finite contribution to the sum over dif- 
ferent sequences of intermediate states a is obtained only 
if the n one-particle hops are such that the final config- 
uration has an overlap with the initial one, correspond- 
ing to the ground state. The n sums over intermediate 
states Qffe in equation ( |l3| ) can then be rewritten as a sum 
over all the sequences S = {ot\, ct2, ■ ■ ■ , Q n -i) which give 
a nonzero contribution. Denoting the numerator of the 



terms by Num(S') and the denominator by Den(S'), equa- 
tion (|13|) takes the form 



E 



(n) 



E 



Num(g) 
Den(5) 



(14) 



We will now evaluate the numerator and the denomina- 
tor separately. The numerator Num(5) can be calculated 
directly from the n hopping matrix elements, thereby tak- 
ing into account the flux dependent phase for hops cross- 
ing the boundary. Since we consider fermions, the corre- 
sponding operators anti-commute and their order in the 
products (H) defining the basis states \ip a ) is crucial for the 
sign. When the starting point (the ground state configu- 
ration) is re-established after n consecutive hopping pro- 
cesses, the order of the operators can be modified. Then, 
the sign of the permutation Ps of the operators, caused by 
the sequence of one-particle hops S 1 , must be incorporated 
in the result. Altogether, one obtains 

Num(S) = sign(P s ) H)" exp [-i<fe(/i £ - &b)] ■ (15) 

hf and hi, denote the number of hoppings across the 
boundary between sites (L x ,y) and (1,2/) in forward and 
backward direction, respectively. Therefore, only the cor- 
rections to the ground state energy due to sequences with 
/if - /ib ^ are flux dependent. 

Moving particles creates defects in the Wigner crystal. 
This increases the interaction energy (see equation |ll|) of 
the ground state E\ by the amount Ue ai where 



(16) 



is non-negative and independent of U . It accounts for the 
difference of the inter-particle distances between the con- 
figurations of the state \ip a ) and the ground state IV'o)- 
We assume always U 3> W, such that the difference in 
potential energy can be neglected in a first step except 
for the case e a = which occurs for some special sample 
geometries and particle numbers. Corr ecti ons due to the 
disorder will be considered in Section (4.2). 

In the generic case when e a > for all intermediate 
states, the energy differences in the denominator Den(S') 
are dominated by the difference in interaction energy. In 
the limit of strong interaction, one therefore gets 



Den(S) » (-UT- 1 J] 



(17) 



a=l 



4 Persistent current 

4.1 Longitudinal current 

Since the longitudinal persistent current at zero tempera- 
ture / = —dEo/d(f> x is given by the flux dependence of the 
ground state energy, a perturbative treatment of the latter 
in terms of the hopping t yields a systematic expansion of 
the persistent current. 
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Fig. 1. 4 particles on 6 x 6 sites. The Wigner crystal con- 
figuration has regular square structure. The arrows indicate 
L x = 6 hops from which a contribution in lowest order to the 
persistent current arises. 



Calculating the n order correction E { 



O), 



to the 



ground state energy, one gets the correction to the persis- 
tent current 



I 



(«), 



dcf) x 



(18) 



in n order in the perturbation H^. 



4.1.1 Relevant terms of the perturbation theory 

When the sequence S of n hopping elements is chosen 
such that each of the particles returns to its initial po- 
sition without completing a tour around the ring, every 
particle which has crossed the boundary must necessarily 
cross it a second time in the opposite direction such that 
hi, = h[ . Therefore the flux dependence of these contribu- 
tions disappears and they cannot influence the persistent 
current. The lowest order of the perturbation theory which 
yields a finite contribution JW is n — L x , corresponding 
to the sequences in which one particle starting at (xo,yo) 
crosses the boundary and returns to its original position 
after completion of its journey around the ring at constant 
y = yo. If there are more than one particle in the line of the 
lattice with constant y = yo, a contribution of the same 
order arises from sequences of hops which move each of 
the particles to the position of its neighbor (see Fig. |l|). 
Since any hopping in y-direction leads to an increase of 
the order of the contribution, the lowest order correction 
to the persistent current is given by the considered pro- 
cesses in which y is kept constant and all of the hops are 
either in forward or in backward direction. 

We now address the dependence of the permutation P$ 
corresponding to this kind of process on the number N y of 
particles in the line at constant y. If N y — 1, the particle 
returns to its initial site, the final configuration is exactly 
equal to the initial one and sign(Ps) = +1. In the general 
case of arbitrary N y , the considered sequences of hops 
(see Fig. [l]) lead to a cyclic rotation of the order of the N y 
particles and the sign of the corresponding permutation is 
sign(Ps) = (-l)^" 1 . 



This yields the result for the lowest-order ^-depen- 
dence of the ground state energy 



E, 



(£»)/ 



E 



^sign(Pg) exp [-i<j) x (hf - h h )} 

U X g £ , , , £ 



Each sequence considered in this sum contains either L x 
forward hops with hf = 1; hb = or it contains L x back- 
ward hops with hi = 0; hb = 1. Each given backward se- 
quence Sb can now be assigned to the forward sequence S{ 
with the reversed order of hops whose contribution differs 
only in the sign of the flux-dependent phase-factor. One 
can express the result as a sum over the forward hopping 
sequences 



E, 



St 



Within the above perturbation theory, when processes cor- 
responding to two loops around the ring (which are at least 
of order 2L X ) are neglected as compared to the lowest or- 
der one-loop processes, the flux dependence of the ground 
state energy is harmonic and 27r-periodic in 4> x . 



4.1.2 Lowest order result for the persistent current 



From (20), one obtains the persistent current in L x -th 
order perturbation theory 



with the flux-independent amplitude 



sign(Ps) 



(21) 



(22) 



Si 



This result contains several interesting features. First, the 
absolute value of the persistent current decays proportion- 
ally to t Lx /C/ Lx_1 , in the limit of strong interaction. In or- 
der to determine the constant prefactor of this power law, 
it is sufficient to figure out all possible processes which 
transform the ground state into itself, using L x forward 
hopping processes, and to calculate the corresponding e a 
from (|lq). 

Furthermore, the sign of the dominating contributions 
to the persistent current in the limit of strong interaction 
is given by — sign(Ps), which for spin-polarized electrons 
is given by (— 1)^" with the number N y of electrons in the 
line of the sample at constant y. This is consistent with 
the well-known theorem by Leggett for the sign of the 
persistent current of spinless fermions in Id [^8| (positive 
/ or paramagnetic response for N even and negative or 
diamagnetic for TV odd). Only if the unperturbed ground 
state (Wigner crystal) configuration of the particles is such 
that the particle numbers N y in different occupied lines y 
have different parity, the prefactors of the corresponding 
terms have to be considered to determine the sign of the 
persistent current. For N y particles in a line of length L Xl 
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Fig. 2. 

lattice. 



Ground state configuration for 4 particles on a 4 x 2 



the number of hopping sequences N seq (N y ) going from the 
ground state to itself is the number of terms contributing 
to the sum in equation (^2|) . Therewith, the result for the 
persistent current can be roughly estimated to be 



(23) 



In the limit of low filling N y /L x — > 0, N seq (N y ) is approx- 
imatively given by the number of possibilities N scq (N y ) 
for N y particles to each make L x /N y (here we assume 
for simplicity that L x /N y is an integer) forward jumps 
to reach the position of its neighbor, leading to the esti- 
mate N scq (N y ) « L x l/[(L x /N y )\] N y. At finite filling, one 
must consider that the neighbor particle must have left its 
starting site before the arriving particle can do its last hop. 
This correlation of the order of the hops of different par- 
ticles reduces N scq (N y ). With increasing filling, N seq (N y ) 
starts to exponentially increase with N y and continues 
to increase more slowly until N y — L x /2 (half filling). 
For larger filling it decreases, thereby obeying a symme- 
try with respect to half filling which is a consequence of 
the symmetry between particles and holes. One can ex- 
pect that the contribution of the line N v — N y lax with 
the largest number A r soq (iV™ ax ) of sequences (at low fill- 
ing this is the one with the maximum number of par- 
ticles) may dominate over the contributions of the ones 
with fewer sequences. The sign of the persistent current is 
then likely to be (—1)^ 



4.1.3 Examples 

As an example, we calculate explicitly the lowest order 
term in jj of the persistent current for 4 spinless fermions 
on a few small L x x L y rectangular lattices, using the 
formula (El. 



4x2 sites We start with the simple case of a 4 x 2 
lattice. The electro-statically lowest energy configuration 
(the Wigner crystal) is shown in Fig. ||. The number of 
particles in the two lines at y = 1 and y = 2 of the sys- 
tem is N y = 2, and the lowest order of the perturbation 
theory which yields a contribution to the persistent cur- 
rent is n = L x = 4. From these two ingredients, we can 
immediately determine the sign (sign(Pg) = — 1 for N y 
even) and the power law of the decrease of the persistent 
current at strong interaction strength. The leading term 



of the amplitude (|22j) is given by 

,4 _ 1 



j(4) 



t 

773 X/ 



Of 



(24) 



This result is always positive and oc l/U 3 . Therefore, the 
response of the system to the applied flux is always para- 
magnetic at strong interaction. 

For this example, it is an easy exercise to figure out all 
hopping sequences which contribute in lowest order. For 
each line y = 1,2, one finds 4 different sequences and can 
explicitly calculate the interaction energies and e Q of the 
intermediate states. This allows to evaluate the sum over 
all sequences in (K4) with the result 



I -» « 853-^- 

U A 



for 1 



0. 



(25) 



In the case of 8 spinless fermions on 4 x 4 sites, the 
same analysis can be carried out except for the fact that 
now four lines y — 1, 2, 3, 4 must be considered, yielding an 
additional factor of 2 in (25). In particular, the persistent 
current is always paramagnetic in the strong interaction 
limit, as found numerically in Ref. [p5[. 



2x4 sites The situation for 4 particles on 2 x 4 sites is 
even simpler since there are four lines y = 1,2,3,4, each 
of them containing N y = 1 particle. Since L x — 2, the 
leading order of the perturbation theory at strong inter- 
action is n = 2, and the sequences of hops contain only 
one intermediate state. The evaluation of (E3) yields 



-15- 



U 



for 



t 

U 



0, 



(26) 



and the persistent current is diamagnetic at strong inter- 
action. 



6x6 sites 4 particles on 6 x 6 sites is the situation in- 
vestigated numerically in Ref. ||, at strong disorder. As 
a function of the interaction strength, an increase of the 
average persistent current at intermediate strength and a 
decrease at strong interaction was found. An exponential 
dependence on the interaction strength was fitted to the 
data for not too strong interaction. Furthermore, it was 
noticed that, at strong interaction, the persistent current 
became paramagnetic for all samples, independent of the 
disorder realization. 

The Wigner crystal ground state on such a lattice is of 
square form, as shown in Fig. |l]. Thus, there are two lines 
with N y = 2 in which hopping sequences of L x = 6 hops 
can contribute to the persistent current. This explains im- 
mediately that the response is paramagnetic, since N y is 
even and — sign(Ps) = +1. The persistent current de- 
creases in the strong interaction limit as t 6 /U 5 . A labori- 
ous evaluation of all the 18 sequences for each of the lines 



with N y = 2 yields the result 



Z-s./W w 1.808 x 10 6 



for 



t 

U 



(27) 
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As compared to the previous cases, the prefactor is much 
larger. This is caused by the bigger number of contribut- 
ing sequences and, more importantly, the difference in in- 
teraction energy between the intermediate states and the 
unperturbed ground state e a . The latter can be very small 
when the distance between the particles is large. 

A numerical investigation by direct diagonalization of 
the corresponding Hamiltonian matrices pQ| confirms the 
signs, the power laws and the numerical prefactors pre- 
dicted by the above formulas, also for the last case of 4 
particles on 6 x 6 sites, where the persistent current is 
indeed found to follow the power law of (^) at strong in- 
teraction. However, it must be noticed that the sign of the 
persistent current is well established Q at interaction val- 
ues much lower than the ones where the agreement in am- 
plitude with our formula starts to be good. Even though 
the data always follow the power laws at strong inter- 
action, fitting an exponential interaction dependence, as 
done in Ref. ||, is possible at moderately strong inter- 
action, and might allow to extract useful informations in 
the regime where higher order terms of the perturbation 
theory are non-negligible. 

4.2 Disorder effects at strong interaction 

The role of the disorder and realization-dependent fluc- 
tuations of the persistent current vanish in the limit of 
strong interaction, when W/U — > 0. In this section, we 
treat the lowest order correction in W/U to the results 
presented above. 

In addition, we consider the special case of a perfectly 
clean system W = 0, in which the translational symmetry 
can considerably influence the interaction dependence of 
the persistent current. 

4.2.1 Disorder corrections to the persistent current 

In order to take into account the disorder corrections in 
the perturbation theory for the longitudinal persistent 
current, it is not sufficient to consider in the denominator 
of the expansion terms the e Q which completely neglect 
the disorder energy. Instead, the full energy differences 
Eq — E a = —Ue a with 

e a = ±{E I a + E»-El-EZ) (28) 

account also for the differences in disorder energy between 
the intermediate states a and the ground state. The dif- 
ference 

?a - e„ = ~ - E») = ^ - > 

(29) 

is of the order W/U and vanishes in the limit W/U — > 0. 
With the definition 

d a = ^ f Vi(n i {a)-n i (0)) (30) 



we can write e a = e a + j^d a and therewith express the 
energy difference terms in the denominator of (|13| ) as 



E -E a Ui a Uea (i + w_ da y 

Taking the lowest order term in W/U, and averaging 
over the ensemble yields 

(snhs)"^ 1 -^ «■>)■ (321 

where the brackets < . . . > denote the ensemble average 
over all disorder realizations. At first glance one could ex- 
pect that the correction linear in W/U vanishes because 
< d a >= 0, when the disorder average is taken over all 
values of d a . However, since the ground state is the Wigner 
crystal pinned at the lowest disorder configuration, d a is 
more likely positive than negative. The first correction is 
thus linear in W, and since e Q is always positive, the lowest 
order correction to the persistent current due to disorder 
is reducing the contributions of all sequences, the result 
being 

~ o fL * sr si s n ( p g) 
w ~ u L *~ l 2^ : ; e 

x(l-E^^^ + 0[(W) 2 ])-(33) 
4.2.2 The particular case W = 

At W — we have Hq = H], and the ground state at 
t = is degenerate because of the translation symmetry. 
The number ud of equivalent degenerate Wigner crystal 
positions and corresponding basis states IV^) depends on 
the system size and the number of particles. The hopping 
terms however lead to a coupling of these degenerate basis 
states, and split the degenerate levels, except for special 
situations where symmetries persist. The coupling terms 
themselves can be expressed in terms of a perturbative 
expansion in t. The order p in which the degeneracy of 
the ground state is lifted may be different from the lowest 
order n which yields a finite persistent current. We address 
in the following the three different cases which may occur. 

p > 7i In this case, the splitting can be ignored as com- 
pared to the persistent current at strong interaction. The 
perturbation theory can be applied as in the disordered 
case, using one of the degenerate ground state configura- 
tions (the result will be the same for arbitrary superposi- 
tions of the degenerate basis states) . 

p = n When the persistent current is given by terms 
which are of the same order as the ones which lift the 
ground state degeneracy, the flux dependent correction to 
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the ground state energy and the persistent current are 
given by the lowest eigenvalue of the effective coupling 
matrix M between the tid degenerate basis states. The 
matrix elements are 



M 0,0' 



E 



(34) 



ai,Q2,---,Ct p _l 



(4 0) \H K \^ ai ){^ ai \H K \iP a2 } . . . (^-Jgxl^ ) 
(E a - E ai )(E — E a2 ) . . . (E a — E ap _ 1 ) 

Since in this case all matrix elements are oc t n /U n ~ 1 , the 
lowest eigenvalue and thus the persistent current will fol- 
low such a power law too. 



p < n When the levels are split at an order which is lower 
than the one which yields contributions to the persistent 
current, the higher order terms must be calculated using 
the lowest energy ground state found from the diagonal- 
ization of M . Such a ground state will be a superposition 



Wi^fptyT) with £lJ>| 2 = 1 (35) 
0=1 

of the different Wigner crystal positions. Plugging this 
ground state into the general expression for the corrections 
to the ground state energy (|l^), one gets 



E^(W = 0) = J2 fpfp 



0,13' 



E 



(36) 



| H K | ) (^ ai | H K | ^ 2 ) . . . (^ a „_ 1 1 H K | ) ) 
(Eq — E ai )(E — E a2 ) . . . (E — E an _ 1 ) 

and realizes that contributions can arise from hopping 
processes starting at any of the ground state components 

\ipo ), and ending at an arbitrary \ijjjf ') where and (3' 
can be equal or different. The processes with f3 ^ f3' can 
in some special cases have a lower power in -j than the 
processes with (3 = (3' , relevant in the disordered case, and 
therefore dominate in the limit of strong interaction. 

The case W = can be qualitatively different from the 
disordered case W ^ 0, even at very weak disorder. The 
reason for this is the fact that the unperturbed ground 
state is a superposition of the degenerate positions of the 
Wigner crystal. This degeneracy is lifted by the disorder, 
and the typical energy difference between two positions is 
of the order Wy/~N (this arises from the N different ran- 
dom on-site energies, all being of order W). The coupling 
between these electro-statically equivalent states is pro- 
vided by the hopping terms and can itself be estimated 
from a perturbation theory. The lowest order coupling 
needs a number of hops which is typically p = N since 
the Wigner crystal must in most cases be translated as a 
whole and each of the particles must hop at least once. In 
some cases, a p < N can be sufficient, as in the example 
of 3 particles on 4 x 2 sites (see below). 

Then, the coupling decreases at strong interaction as 
t p /J7 P_1 and becomes, in the limit of strong interaction 



Fig. 3. Two degenerate configurations for 3 particles on 4 x 2 
sites which can be linked to each other by sequences of only 
two hopping processes. 



U, always much smaller than the splitting due to the dis- 
order. This holds at any finite disorder value (provided 
p > 1 ) , and prevents a mixing of the different Wigner crys- 
tal configurations in the ground state. Only for W = 0, 
the splitting by the disorder vanishes, and an infinitesi- 
mal coupling leads to the superposition of the different 
configurations. 



4.2.3 Example for disorder corrections 

4x2 sites, 4 particles The disorder dependence of the 
flux sensitivity is illustrated for the example of 4 particles 
in (L x ,L y ) — (4,2). One expects an 1/U 3 dependence for 
the lea ding term of the persistent current (pj) and, accord- 
ing to (p3[), a W/U 4 dependence for the linear correction 
in disorder. We have performed numerical calculations at 
strong interaction |40| , which are in good agreement with 
these power laws predicted by the lowest order of the per- 
turbation theory. The persistent current decreases when 
the strength of the disorder is increased, as predicted by 
the negative correction in formula (|33"|). 



Effects of degeneracies at W = The clean case how- 
ever has a very different amplitude due to the degeneracies 
in the Wigner crystal. The persistent current can be calcu- 
lated from ([56]), and after considering all the possible se- 
quences between the different degenerate configurations, 
one obtains I^(W = 0) w 2.1238 j( 4 \ with i"( 4 ) taken 
from (^4|). A numerical check |40 confirms this result at 
strong interaction. 



4x2 sites, 3 particles, W = In order to show that 
the absence of disorder does not only change the prefac- 
tor of the interaction dependence at large U, but can even 
influence the power of the [/-dependence, we consider 3 
particles in (L x ,L y ) — (4,2). The particle which is alone 
in one line has two possible positions which are degenerate 
at t = W = 0, and there exists a sequence of only two one- 
particle hops which links the two degenerate states with 
each other (see Fig. ^). Since the ground state at W — is 
a superposition of all possible degenerate configurations, 
this process leads to a contribution oc t 2 /U to the persis- 
tent current which dominates at strong interaction. In the 
disordered case, in contrast, the degeneracy of the ground 
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state configurations is lifted and only sequences consist- 
ing of at least four hops can link the ground state to it- 
self. Then, the interaction dependence of the lowest order 
contribution to the persistent current is oc t A / (WU 2 ). In- 
terestingly, in this example, one of the intermediate states 
touched by a sequence of 4 hops can be a different, but 
electro-statically equivalent configuration of the Wigner 
crystal, and e a = 0, such that even at U ^> W, the disor- 
der energy cannot be neglected. 

A numerical check pQ | shows that by adding the dis- 
order, the asymptotic dependence indeed follows this pre- 
diction and goes from oc U" 1 to oc U~ 2 . 



4.3 Transverse current 

4.3.1 General considerations 



Fig. 4. A sequence of hopping processes which gives rise to a 
finite contribution to the transverse current in the disordered 
system of 4 particles on 6 x 6 sites. The arrows indicate L x + 
L y — 12 forward and upward hops. 



Without disorder, the current in y-direction must vanish 
because of the symmetry of the problem (no flux is applied 
in y-direction) . While this remains true in the ensemble 
average, the presence of disorder breaks this symmetry 
in individual samples. A finite transverse current can be 
observed in addition to the longitudinal current driven by 
the flux <j) x . The transverse current can be calculated from 



dE ( 



(37) 



„=o 



L y ('down- 



such that a theoretical approach to determine It(4>x) con- 
sists in using the perturbative expansion for the ground 
state energy, as for the longitudinal current, but with a fi- 
nite 4> y . Then, the derivative with respect to 4> y at <p y = 
can be evaluated for the leading contribution in the strong 
interaction regime. Therefore, only terms in the series (O) 
which depend on (j> y can make a finite contribution to the 
transverse current. 

In order to get such contributions, one needs hopping 
sequences in ^-direction which contain hops crossing the 
boundary between sites (x, L y ) and [x, 1) with h u —hd =/= 0. 
h u and /id denote the number of such hops from y = L y 
to y = 1 ('upwards') and from y = 1 to y 
wards'). 

However, the flux dependence of the energy due to se- 
quences containing L y upward or L y downward hoppings 
is oc cos 4> y , as for the longitudinal processes (^p|), which 
are oc cos <p x . As a consequence, the resulting transverse 
current from these sequences is oc sin (f> y , and vanishes 
when <p y = is taken, just like the longitudinal current 
( ^l| ) at (j> x = 0. Therefore, no transverse current can be 
obtained in the order L y of the perturbation theory. 

Nevertheless, higher order sequences exist which do 
lead to finite contributions. These are due to sequences of 
hops which cross both, the boundary in x, and the bound- 
ary in y-direction. The lowest order of the expansion in 
which transverse currents appear is to = L x + L y , cor- 
responding to sequences which consist of a combination 
of only forward and downward hops, only backward and 
upward hops, or vice versa, which complete one round in 



both spatial dimensions of the system. An example of such 
a sequence is shown in Fig. |J 

As before, each of the sequences containing only for- 
ward and upward hops 5f u can be associated with the re- 
verse process consisting of backward and downward hops 
Sbd- The contribution of the reverse process differs only 
in the sign of the flux-dependent phase factors. The cor- 
rection to the ground state energy ( [l3"| ) of all Sf u and Sbd 
sequences together can be written as 



E, 



(m) 

0,fu+bd 



u r - 



— COS((f> X +(f)y) 



sign(P gf J 



(38) 



The same association can be made for sequences contain- 
ing only forward and downward hops S^, which are re- 
grouped with the corresponding reverse processes consist- 
ing of backward and upward hops, yielding 



p( m ) _ _o 
-^O.fd+bu ~~ * jj 



jCOS^-^y)^ _TSrr_g 
C ai . . . 6 Q 

i>fd 



;D(PsJ . (39) 



4.3.2 Lowest order result for the transverse current 

By taking the derivative with respect to <j) y at (f> y 
we obtain for the transverse current in TO-th order of the 
perturbation theory 



= 0, 




sm < 



(40) 



sign(Ps fd ) 



with m = L x + L y being the lowest order which yields a 
non-zero contribution. Even when U 3> W, the disorder 
energy cannot be neglected as in the lowest order term 
for the longitudinal current, since the two contributions 
to the current in (^) cancel each other exactly when the 
disorder vanishes. It can also be seen that the lowest order 
result for the transverse current is proportional to sin^ 
and vanishes with the longitudinal current at <p x = 0. 
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4.3.3 Example 

In the case of 4 particles in a lattice of 6 x 6 sites, we 
have seen above (see equation (§7|)) that the longitudinal 
persistent current decreases at strong interaction following 
the power law oct 6 /U 5 . The suppression of the transverse 
current is much more pronounced, since the leading contri- 
bution at strong interaction is of the order L x + L y = 12, 
and our theory ( |40| ) predicts that it decays oc t 12 /!/ 11 . 
This means that the particle mobility is more and more 
restricted to the longitudinal direction when the interac- 
tion is increased. That the suppression of the transverse 
current is much stronger than the one of the longitudinal 
persistent current has been noticed in the numerical study 
of Ref . H . The dominance of the longitudinal current at 
strong interaction also explains the observation of "plas- 
tic flow" in a study of the local persistent currents |$5|. 
The orientation of the local currents has been proposed as 
a signature of the insulator-metal transition occurring in 
interacting 2d systems || . 



4.4 Persistent current for electrons with spin 

In this section we present some examples of the effects 
of the spin of the electrons on the longitudinal persistent 
current, showing that our approach is not restricted to 
spinlcss fcrmions, but can also be used for electrons car- 
rying spin. 

For finite systems containing N particles, it is well 
known that the spin polarization of the ground state can 
depend on the magnetic flux due to level crossings between 
states of different spin symmetry as a function of the flux 
. However, without an external magnetic field acting 
on the electrons in the ring, there is a degeneracy related 
to the operator S z in the subspace of fixed total spin and 
one can write E(S 2 , S z ) = E(S 2 ). 

It is in principle possible (though difficult experimen- 
tally) to create a magnetic flux through the ring while the 
magnetic field remains vanishing at the positions of the 
electrons in the ring. Such a flux does not lift the spin 
degeneracy mentioned above. If we now want to follow 
the dependence of the ground state energy on the flux, we 
can restrict the study to the subspace of minimum abso- 
lute value of S z , choosing S z = for an even number of 
particles with JV-f = N/2 spins up and N± — N/2 spins 
down, or S z — 1/2 for an odd number of particles with 
iVj = spins up and N[ = spins down. 



4.4.1 One-dimensional systems 

In Id, a general rule exists for the sign of the persistent 
current in the case of spinless fermions p£| ]. This rule is 
valid at arbitrary disorder and interaction. Below, we ad- 
dress the question of the sign of the persistent current for 
electrons with spin in different parameter regimes. 



Non-interacting electrons with spin In the case of non- 
interacting electrons with spin, one can separately con- 
sider the flux dependence of the energy of the particles 
with up spin and the one of the particles with down spin, 
and add up their contributions to the persistent current. 

If N is even, the contribution of the N/2 up-spins to 
the persistent current will have the same sign (— 1) N / 2 as 
the contribution of the N/2 down-spins. This determines 
the sign of the persistent current. For odd N, the con- 
tribution of the (N + l)/2 up spins will have the sign 
(— 1)( JV + 1 )/ 2 while the down spins contribute a term with 
sign (— l)^ -1 )/ 2 . In this case, to know the sign of the per- 
sistent current one must compare the amplitudes of the 
two contributions and the result will in general depend on 
the disorder configuration. However, in the zero-disorder 
case at low filling, the sign of the persistent current around 
flux (f> = is known to be paramagnetic for odd N > 1 

& 

iV even, strong interaction At strong interaction, the 
charges form a Wigner crystal which is pinned by the dis- 
order. In this limit, the spin degree of freedom can be 
treated by an effective spin Hamiltonian. The latter turns 
out to be the anti-ferromagnetic Heisenberg Hamiltonian 
when = 0, the positions i = i±, 12, ■ ■ ■ , iff of the charges 
of the Wigner crystal being the spin lattice sites. Thus, 
the expectation values of the occupation number operators 
vanish everywhere except on these sites where they satisfy 
< ^»*,T"I"^»*4 Fb r a- n even number N of spins in Id, 

according to a theorem by Marshall, the ground state |^o) 
of this spin Hamiltonian is a singlet of total spin S = 
Q , and can be expressed in the form 

i*>=EAn(<t)"" ,m n(< i )""' m io) <«> 

fj i—ii i=ii 

with real fg > for all spin configurations (3 with fixed 
S z = 0. Note that in this expression, the ordering of the 
operators is done firstly according to spin, and secondly 
according to the position. 

As in the case of spinless fermions without disorder, 
the ground state is a superposition of different basis states 
and the lowest order contributions to the ground state 
energy, which are flux dependent, are again of the order 
n = L x . Similar to (|36|), they can arise from different 
kinds of sequences. First, only the up spins are moved 
around the ring, giving rise to a cyclic permutation of the 
spin up operators in (f4l|), yielding the sign (— 1)^. The 
sequences involving only down spins give the same result. 
In addition, a sequence which moves all of the particles to 
the position of their neighbor can also contribute to the 
flux dependence since the resulting spin configuration is 
also contained in the ground state. The sign however is 
the same as the one of the previous sequences, since only 
one of the particles crosses the boundary and therefore 
only the order of the operators for one spin direction has 
to be restored by the corresponding cyclic permutation, 
and because the prefactor fgfgi, arising from the weights 
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of the different components in the ground state, is always 
real and positive. 

Since all contributions have the same sign, the sign 
of the persistent current in one-dimensional disordered 
chains at strong interaction is always given by this sign 
(—1) N / 2 , provided the particle number is even, indepen- 
dent of the particle density N/L. Because the ground state 
( fill ) holds only at <fi = 0, this sign rule is granted only in 
the vicinity of 4> = 0. 

The resulting sign is the same as the one found in the 
non- interacting case [Q, and for the Hubbard model at 
half filling |3l| , but it differs from the result for the Hub- 
bard model at low filling (3^,^], where the current around 
4> = is found to be diamagnetic for even numbers of par- 
ticles at strong interaction. This difference may be due to 
the fact that the Hubbard interaction is local, making the 
Hubbard model at low filling in the U — > oo limit equiv- 
alent to spinlcss fermions. The more realistic long-range 
Coulomb interactions considered here do not show this ar- 
tifact and lead to a qualitatively different behavior. The 
result is also different from the one obtained for Coulomb 
interacting electrons with spin in a continuous model with 
one barrier [B3J , which is always found to be diamagnetic 
at strong interaction. In this case, one may attribute the 
result to the fact that the interaction could lead to a rigid 
Wigner crystal which, in a continuous model, would be 
equivalent to one single heavy particle. 



4.4.2 Two-dimensional systems 

For two dimensional systems, there is the same degeneracy 
in the position of the up and the down spins in the mini- 
mum S z subspace. In the strong interaction limit, each line 
of the Wigner crystal should follow the one dimensional 
law explained above, but the different spin configurations 
contributing to the ground state in 2d may have differ- 
ent numbers of spin up and spin down electrons in a given 
line. This yields contributions of different signs to the per- 
sistent current which have to be considered explicitly. It 
seems not to be possible to provide a simple sign rule. 
However, the power law decay of the persistent current at 
large interaction strength should have the same exponent 
as in the spinless case. 



4.5 Size dependence and limitation of the theory 

Here, we address the size dependence of the persistent 
current and the limitation of our theory which is given by 
the presence of defects in the Wigner crystal. 



4.5.1 Localization length for the Wigner crystal 

From the exponential size-dependence of equation (p2|), 
one can extract the localization length for the Wigner 
crystal. A comparison to the exponential size scaling 



yields for large L x a localization length 



(43) 



showing the same decrease with increasing interaction 
strength as the localization length of the Mott insulator 
appearing in the Hubbard model at strong interaction p|| . 



4.5.2 Defects in the Wigner crystal 

In the thermodynamic limit, the Wigner crystal is no 
longer a single domain and the above theory cannot be 
applied in the present form. At any finite, even very strong 
interaction, defects and domain walls can arise which al- 
low to gain an amount of potential energy which grows 
with the size of the domains. The crystal might then pre- 
fer to be divided in defect-free domains each being pinned 
by the disorder. In the weak disorder limit, one can use 
similar arguments as in in order to estimate the size 
of such domains. In the following, we briefly address dif- 
ferent types of defects and calculate the associated critical 
ratio (jjp) c above which the theory applies. 



Point defects In the case of point-like defects, consisting 
of a single charge of the Wigner crystal being deplaced 
by one lattice constant a with respect to its position in 
the perfect crystal, the gain in disorder potential energy 
can be estimated to be AEjj°™ der ~ W. On the other 

hand, the cost in interaction energy is Afteraction ~ 
with b being the mean spacing between particles. Since 
the charge density is given by v = 1/b 2 , this yields the 
density depending criterion 



U_ fU_ 

w > V w 



point 



l 



(44) 



for the stability of the Wigner crystal against the creation 
of point defects. 



Domain walls The creation of a domain wall costs an 
interaction energy which is of the order Zi£'™g raction w 

^L wa11 with L wa11 being the length of the wall. By de- 
placing a domain of linear size R containing Nd = R 2 v 
particles, one can typically gain an amount W\fWd~ of dis- 
order potential energy. However, the first domain-like de- 
fect will appear in the most favorable position. Since the 
maximally possible gain is the extremely unprobable value 
WNd, one can speculate that the gain in the optimal posi- 
tion is given by an intermediate power AE^?™^ ss WNj 
with 1/2 < 7 < 1. The length of the corresponding wall 
around the domain is £ wa11 ~ \/N2 such that this results 
in the stability of the crystal against the creation of do- 
mains of size R for 



I(L X ) cx exp(-L x /£) 



(42) 



U 

TP 



> 



domain 



(R) 



R 2f- 



(45) 
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At a fixed W, if U > U* omaia (L), a crystal of size L will 
not be perturbed by domains. Otherwise, one can use U = 
C/ c domain (i? c ) to extract the typical domain size R c . 

In systems larger than R c , the electron crystal is di- 
vided in many domains of size R c . If we neglect the cou- 
pling energy between these domains, a rough estimate for 
the response to a flux <f> x in ^-direction of a system of size 
{L x ,L y ) is to take the product of the amplitudes of the 
longitudinally aligned L x / R c domains and to sum over the 
responses of the L y /R c "channels". 

Although we have not addressed all possible kinds of 
defects, we can assume that there is always a critical 

threshold (^) above which our perturbative theory 
applies. 



5 Summary 

We have presented a systematic perturbative treatment 
of persistent currents in 2d lattice models for the case 
of strong Coulomb interaction, when the electronic charge 
density forms a Wigner crystal. The contribution with the 
weakest interaction dependence corresponds to sequences 
of one-particle hops along the shortest paths around the 
ring. These sequences dominate in the limit of strong in- 
teraction and determine the sign of the persistent current. 
For spinless fermions, this sign follows simple rules which 
depend only on the structure of the Wigner crystal. 

Furthermore, we have shown that, except special cases, 
the leading term for the persistent current and therewith 
the sign of the persistent current at strong interaction do 
not depend on the realization of the disordered potential. 
Only the complete absence of disorder can qualitatively 
change the behavior. We considered the disorder correc- 
tions systematically and showed that they decay as W/U 
at strong interaction. 

In addition, we have shown that transverse currents 
appearing in individual disordered samples are suppressed 
much faster than the longitudinal current, thereby estab- 
lishing an orientation of the local currents in longitudinal 
direction when the interaction is increased. This explains 
the numerical observation of the realization-independent 
sign of the persistent current and of the orientation of the 
local currents in longitudinal direction reported in Refs. 

linn. 
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